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Abstract. There are indications that interstellar and interplanetary dust grains have an inhomogeneous and fluffy 
structure. We investigate different methods to describe light scattering by such composite particles. Both a model 
of layered particles and discrete dipole calculations for particles with Rayleigh and non-Rayleigh inclusions are 
used. 

^ . The calculations demonstrate that porosity is a key parameter for determining light scattering. We find that the 

' optical properties of the layered particles depend on the number and position of layers if the number of layers 

l/^ , is small 15). For a larger number of layers the scattering characteristics become independent of the layer 

■ sequence. The optical properties of particles with inclusions depend on the size of inclusions provided the porosity 
0^ ' is large. The scattering characteristics of very porous particles with inclusions of different sizes are found to be 

, close to those of multi-layered spheres. 

' We compare the results of these calculations with the predictions of the effective medium theories (EMT) which 

are often used in astronomy as a tool to calculate the optical properties of composite particles. The results of 
. our analysis show that the internal structure of grains (layers versus inclusions) only slightly affects the optics of 

CL(' particles provided the porosity does not exceed 50%. It is also demonstrated that in this case the optical properties 

JL I of composite grains calculated with EMT agree with the results of the exact method for layered particles. For 

■ larger porosity, the standard EMT rules (i.e., Garnett and Bruggeman rules) give reliable results for particles 
' with Rayleigh inclusions only. 
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H 1- Introduction Fortunately, scientists had felt the necessity to treat 

- - ' the scattering by composite and inhomogeneous particles 

Various processes operating in interstellar and circum- ^^^-^ consisting of several components even earlier 

stellar media are believed to produce mhomogeneou^^^ ^^^^ existence of interstellar dust was established, 

porous cosmic dust grains (DorschnerCMi Drained . q^^^^^^ (im was the first to find the averaged (effective) 

However, the real structure of interstellar grains remains dielectric functions of such a medium assuming that one 

to be established. The particles may have a layered struc- ^^^^^.^^ ^ ^^^^.^ ^j^^^^ material) in which another 



ture because of the their formation in circumstellar en- 



material was embedded in the form of small inclusions 



vironments (Dommik et al. and further evolution ^^u^j Maxwell-Garnett mixing rule of the Effective 

in molecular clouds (Ehrenfreund[1999, Greenberg|T999). y^^^^^^. Theory; EMT). Later, Bruggeman mB de- 

On the other hand, collisions of grains tend to induce co- ^^^^^ ^^^^^^^ ^^^^ ^^.^-^ symmetric with respect to 

agulation and partial destruction of particles. Because of ^-^^ niaterials. These classical mixing rules are still the 

this, interstellar grains should have inclusions of different ^Qgf^ popular ones 
size. Therefore, considering the optical properties of cos- 
mic dust grains we are forced to solve the difficult problem 

of the interaction of radiation with composite particles of 

j.rr , , , good accuracy, 

ditterent structure. ° 

This first became possible at the beginning of the 

Send offprint requests to: N.V. Voshchinnikov, nvv@astro. 1950s when the Mie solution for homogeneous spheres 

spbu.ru was generalized to core-mantle spherical particles in 



Many scientific and applied problems require calcula- 
tions of light scattering by inhomogeneous particles with 
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three independent papers (Aden & Kerker [T951I Shifrin 
HM^ Giittler IDB^ . Guttler's solution was used by 
Wickramasinghe (1963) who first calculated the extinction 
of layered (graphite core-ice mantle) analogues of cosmic 
grains. 

Mathis & Whiffen lfTU5g|l introduced the first consis- 
tent model of composite cosmic grains which were very 
porous (the volume fraction of vacuum 80 %) aggregates 
of small amorphous carbon, silicate and iron subparticles. 
The optical properties of such particles were calculated 
with the Mie theory and EMT. Mathis l(TM^ updated 
the composite grain model taking into account the abun- 
dances of heavy elements obtained for cluster and field 
B stars and young F, G stars (Snow & Witt 1996). The 
new model consisted of three dust grain populations where 
the visual/near-IR extinction was explained by aggregates 
with ^ A5% vacuum in volume. 

Now light scattering computations for inhomogeneous 
(composite) particles with layers or inclusions from differ- 
ent materials or aggregate particles are often made using 
the discrete dipole approximation (DDA) or simpler the- 
ories (see Voshchinnikov .2004i for discussion) . Note that 
the DDA is a method which still is computationally de- 
manding. Therefore, it is mostly used for illustrative cal- 
culations and not for mass production (e.g., Wolff et al. 
[BMIITMHI Vaidya et aLEOHIl Andersen et al. nUU^ . 

The idea to represent composite interstellar grains by 
multi-layered spheres (Voshchinnikov & Mathis [19991 see 
also lati et al. I2001|l has no immediate physical justifica- 
tion (although such particles may form in stellar envelopes 
and molecular clouds), but is very attractive as an exact 
method to calculate the optical properties of composite 
particles. Such a model permits us to include an arbitrary 
fraction of any material, and computations do not require 
large resources. However, the distribution of material in- 
side particles is always spherically symmetric even when 
its volume fraction tends to zero. 

In this paper, we compare the optical properties of 
composite interstellar grains of various porosity obtained 
from calculations for layered spheres, pseudospherical par- 
ticles with inclusions and homogeneous spheres with an 
effective refractive index. The description of the parti- 
cle models is given in Sect. |21 We compute different ef- 
ficiency factors, albedo, etc. and analyze how these quan- 
tities depend on the order and number of layers and the 
size of inclusions (Sect.OJ. Special attention is paid to the 
consideration of very porous grains fSect. because of 
their particular importance in astronomy, for example, for 
modelling of comets (Greenberg & Hage I1991j) and the 
disc of (3 Pictoris (Li & Greenberg I1998|l . The possibility 
to describe the light scattering by porous particles using 
Mie theory with different EMT rules is studied in Sect.0] 
Concluding remarks are presented in Sect. El 

2. Models of composite grains 

Processes operating in the winds of late-type stars such as 
grain nucleation and growth, shattering in the diffuse in- 



terstellar medium, and finally coagulation and accretion in 
molecular clouds and protoplanetary disks lead certainly 
to dust grains with rather irregular shapes and very com- 
plicated internal structure (Dorschner & Henning I1995|l . 
The details of the grain interiors are not directly impor- 
tant for surface chemistry, but the optical behaviour of the 
particles may be a strong function of this structure. Direct 
evidence for the structure of these particles is difficult to 
obtain with the exception of interplanetary dust grains 
collected in the upper atmosphere and the solar system. 
Therefore, a more general attempt to describe the parti- 
cles and to explore changes in their optical properties is 
required. 

A frequently used approach in astronomy is the mod- 
elling of inhomogeneous grains by two-layered (core- 
mantle) spheres and particles with voids or inclusions us- 
ing EMT-Mie calculations. In this paper, we consider lay- 
ered particles and particles with inclusions as the models 
for the description of the optical properties of the inho- 
mogeneous or composite grains. The amount of a material 
in such particles is determined by its volume fraction Vi 
(S,F,/Vtotai = !)■ The particle porosity P (0 < P < 1) is 
introduced as 

r = Kac/14otal = 1 - 't^solid/T4otal, (1) 

where Vvac and 14oiid is the volume fraction of vacuum and 
solid material, respectively. 

In our calculations presented below, composite par- 
ticles of several materials are considered. The refractive 
indices for them were taken from the Jena-Petersburg 
Database of Optical Constants (JPDOC) which was de- 
scribed by Henning et al. H1999|l and Jager et al. H2UU3|I . 

Carbon and silicates are the materials most often 
used in cosmic dust models (see Mathis et al. 119771 
Draine & Lee 1198^1 and so on). We consider the parti- 
cles composed of amorphous carbon (ACl), astronomical 
silicate (astrosil) and vacuum with varied volume frac- 
tion of each constituent. The optical constants for ACl 
(to = 1.98 -I- 0.23i) and astrosil (to = 1.68 -I- O.OSi), cor- 
responding to the wavelength A = 0.55 /im, were taken 
from the papers of Rouleau & Martin H1991|l and Laor & 
Draine H1993|l . respectively. 

2.1. Layered particles 

Here we represent composite grains by particles consisting 
of many concentric spherical homogeneous layers of cycli- 
cally changing materials. Such a model is not primarily 
meant as a physical description of the actual grain struc- 
ture, but as a possibility for describing light scattering by 
inhomogeneous particles of complex structure. 

Vacuum can be one of the materials, and a composite 
particle may have a central cavity or voids in the form of 
concentric layers. This model allows one to include at any 
position inside a spherical particle any fraction of a mate- 
rial (from extremely small to very large). The light scat- 
tering calculations are based on the exact theory which is 
true for particles of any size and refractive index. 
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Fig. 1. The cut of the spherical particles by the plane. 
The models of composite grains containing the same 
amount of carbon and silicate are shown. The 3- and 9- 
layered spheres consist of equivolume spherical layers with 
the total volume fractions of carbon, silicate and vacuum 
equal to 33.33%. The core-mantle particle includes the 
same mass of carbon and silicate but is free of vacuum. 



The schematic representation (cross cut) of layered 
spherical grains is given in Fig. ^ The particles are com- 
posed of a specified number of concentric spherical homo- 
geneous layers. The order of the layers and their total num- 
ber can be specified separately. Following Voshchinnikov 
& Mathis we assume further that different mate- 

rial layers cyclically change inside a particle and call the 
repeating set of such layers a shell. This means that a par- 
ticle consists of a specified number of concentric shells, and 
the simplest model particle contains one shell of two ma- 
terials. The core-mantle particle presented in Fig.^ does 
not contain vacuum, but its mass is the same as that of 
the other two particles shown. As a result, its volume is 
less by factor of 1/3 and the outer radius by ^1/3 ~ 0.69, 
respectively. 

The formal solution to the light scattering problem for 
n-layered spheres can be easily written in matrix form 
with the separation of variables method (see, e.g., Kerker 
I1969|l . In this case the scattering coefficients are calculated 
as the ratios of two determinants of order 2n+l containing 
Riccati-Bessel functions and their first derivatives of real 
and complex arguments. However, for practical reasons, 
it is better to use the recursive algorithm developed by 
Wu & Wang 1)1991(1 and Johnson H1996(l . In order to make 
calculations for highly absorbing particles of large sizes, 
one should take into account the modifications suggested 
by Wu et al. lfTM7)l and Gurwich et al. ((21111) ■ 

2.2. Particles with inclusions 

The optical properties of particles with inclusions can be 
estimated on the basis of rather complicated calculations 
(see, e.g., Wolff et al. ITTMI Videen & ChylekCMHI) or 
from laboratory measurements (Kolokolova & Gustafson 
121)01) . If the volume fraction of inclusions is not very 
large 10%), the EMT-Mie calculations give the results 
with good accuracy (Wolff et al. IBMl ITMHl Kolokolova 
& Gustafson .200 1,1 . 

However, our goal is the consideration of particles with 
an arbitrary amount of inclusions and different porosity. 
Therefore, the calculations are performed with the discrete 



dipole approximation (DDA). We use the last version of 
the DDA program (DDSCAT 6.0) developed by Draine 
& Flatau (2003). This technique can treat particles of ar- 
bitrary shapes and/or of inhomogeneous structure. A de- 
tailed review of the DDA and its applications is given by 
Draine iPTOjl . 

The particles ("targets" in the DDSCAT terminology) 
are constructed using two special routines. One routine 
produces spherical targets with inclusions of a fixed size, 
while another creates targets with a given distribution of 
inclusions over their sizes. Both routines produce first a 
cube with randomly distributed cubic inclusions. The sizes 
of the target dmax and of the inclusions c?inci are expressed 
in units of the interdipole distance d. In the cube, we just 
inscribe a sphere and remove all inclusions and their parts 
being outside the sphere. 

In contrast to previous modelling (Henning & 
StognienkoCnni Lumme & Rahola 1994, Wolff et al.HnSl 
I1998L Vaidya et al. 12001(1 . porous particles are not pro- 
duced by removing dipoles or inclusions from a target but 
by attributing the refractive index m = 1.000001 -t- Q.Oi to 
them. We believe that such a structure better corresponds 
to cosmic aggregates. 

For the purpose of treating very porous particles, the 
number of dipoles in pseudospheres is taken quite large. 
In all cases considered, the particles with the maximum 
size o?max — 91 are studied. This value corresponds to 
the total number of dipoles in pseudospheres A^dip = 
357128 — 381915 depending on the size of inclusions dinci- 
Thus, the criterion of the validity of the DDA for ex- 
tinction/scattering cross sections |TO|/c(i < 1^ of Draine 
& Flatau ( 2TO3) (fc = 27r/A is the wavenumber with A be- 
ing the wavelength in vacuum) is satisfied up to the size 
parameter soporous = 27rr/A « 27 — 28. 




Fig. 2. The cut of the pseudospherical particles with the 
maximum size c?max = 45 by the plane. The models of 
very porous particles with small single size (left) and dif- 
ferent size (right) cubic inclusions are shown. The volume 
fractions of carbon and silicate are equal to 5%. 



Note that sometimes the fulfillment of this criterion does 
not guarantee the correct result (see, e.g., discussion in 
Andersen et a.l. TM^ . 
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Targets with the values of dinci ranging from 1 to 9 are 
considered. The resulting structure (cross cut) of pseudo- 
spherical grains is shown in Fig. [3 Note that the inclusions 
of the size di^ci = 1 are dipoles, while the inclusions with 
t^inci = 3, 5, 7 and 9 consist of 27, 125, 343 and 729 dipoles, 
respectively. 

The optical characteristics of pseudospheres with in- 
clusions were averaged over three targets obtained for dif- 
ferent random number sets. The calculations show that 
in our case such an approach is practically equivalent to 
time-consuming numerical averaging over target orienta- 
tions. 



3. Towards unified optical properties 

It is evident that composite particles of various structure 
and shape should exist in space. Here we focus our at- 
tention on spheres of different porosity because of their 
frequent use in the modelling of interstellar and cometary 
grains and possible importance in attacking the problem 
of cosmic abundances. Our analysis should also help to es- 
timate the range of validity of some previous models and 
to clarify the physical background of them. 



3.1. Particles of moderate porosity 

We start with the consideration of layered particles. 
Figure |3I shows the extinction efficiency factors Qcxt — 
Cext/T^fg (C'oxt is the extinction cross section, the outer 
radius) of layered spheres. The optical properties of core- 
mantle spheres have been studied rather well and seem to 
show no significant peculiarities^ (Babenko et al. I2003II . 
In contrast, three-layered spheres already can produce 
anomalous extinction of light. The order of the materi- 
als strongly affects the behaviour of extinction for such 
particles (the upper panel of Fig. ^ . First of all, the loca- 
tion of vacuum (the core or the middle layer) is important. 
The curve for particles with a carbon core and an outer- 
most astrosil layer is the most peculiar curve. Here, a very 
rare situation is observed: the first maximum is damped, 
but there is a very broad second maximum. Note that the 
scattering efficiency depends more strong on the order of 
layers than the absorption efficiency. However, all the pe- 
culiarities disappear when the number of layers increases: 
the difference between the curves becomes rather small for 
particles with 9 layers (3 shells) and is hardly present for 
particles with 18 layers (6 shells; see Figs. O and QJ. Figure 
0] shows the size dependence of the scattering (Qsca) and 
absorption (Qabs) efficiency factors, albedo A = Qsca/Qcxt 



Except for a resonance peak arising for particles with 
mantles having a large refractive index (see, e.g., Gurwich et 
al. .2U03J . 
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Fig. 3. Size dependence of the extinction efficiency factors 
for layered spherical particles. Each particle contains an 
equal fraction of amorphous carbon (ACl), astrosil and 
vacuum (the porosity V = 1/3) separated in equivolume 
layers. The cyclic order of the different material layers 
is indicated (starting from the core). The effect of the 
increase of the number of layers is illustrated. The thick 
line at the lowest panel corresponds to compact spheres 
consisting of ACl and astrosil. For a given value of the 
size parameter, the compact and porous particles have the 
same mass. 



and the parameter g^^"^ describing the asymmetry of the 
phase function F{Q, $) 
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(cose) 



/^^f(e,$) cosedc^ 
/4^F(e,<i>)dc. 



(2) 



for multi-layered spheres. As has been noted by 
Voshchinnikov & Mathis H1999fl : the optical properties 
weakly depend on the order of materials in each of the 



The notation g is related to the Henyey-Greenstein phase 
function that is very often use in radiative transfer modelling. 
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shells and are close to some "average" properties, if a par- 
ticle consists of many shells (^3 — 5). In other words, for 
such particles, different efficiency factors as well as albedo 
A and the asymmetry parameter g depend practically only 
on the volume fraction of materials. 

The solid thick lines in Figs. |21 and ^ show the size 
dependence of the optical characteristics of compact {V = 
0) spheres consisting of the same amount of solid materials 
as the porous ones. To compare the optical properties of 
porous and compact particles, it is helpful to normalize 
the size parameter of either the compact or porous particle 
using the relation 

-^compact ^compact /q\ 

In the case of the particles presented in Figs.|31and01 this 
leads to stretching of the x scale for compact particles by 
a factor of •y/3/2 « 1.145. It can be seen that the presence 
of vacuum inside the composite particles reduces the peak 
of the absorption efficiency (Fig. 0] left lower panel) and 
shifts that of the scattering efhciency (Fig. ^ left upper 
panel). Correspondingly, these two effects explain the be- 
haviour of the curves for the extinction (Fig.O. A medium 
porosity influences the albedo and the asymmetry param- 
eter only in a restricted range of size parameters. 

Similar calculations for particles with inclusions were 
performed many times using the DDA technique. But so 
far only Wolff et al. (1994, 1998) considered particles with 
large (non-Rayleigh) inclusions. They computed the opti- 
cal properties of silicate spheres and spheroids with a size 
parameter up to Sporous = 10 and a volume fraction of 
vacuum inclusions up to 80%. For the efficiency factors 
and asymmetry parameter, the difference between parti- 
cles with Rayleigh** and non-Rayleigh vacuum inclusions 
became noticeable for the porosity V ^ 0.4. 

The results of our DDA calculations of the extinc- 
tion efficiency factors for pseudospheres with the porosity 
V = 0.33 are shown in Figs. EHZI The volume fractions of 
the materials are approximately the same as in the par- 
ticles presented in Figs. 13 and ^ but the materials are 
present in the form of inclusions instead of layers (vac- 
uum is considered as a matrix). Figure shows the results 
obtained for particles with single size inclusions. Note that 
despite the different structure of the targets (the number 
of inclusions reduces from ~ 240000 to ~ 330 with grow- 
ing ndip), the differences between the extinction efficien- 
cies are quite small. The same conclusion is correct for 
particles having inclusions of different sizes. The results 
for three different targets with a distribution of inclusions 
of five sizes are shown in Fig. El The size of the inclusions 
ranges from 1 to 9 and the number of inclusions is in- 
versely proportional to their volume. In other words, the 
total numbers of dipoles in the inclusions of each size are 
approximately the same. 

The extinction efficiencies of the particles with inclu- 
sions are compared with those of layered particles in Fig.[71 

* In this case, the size of the inclusion is much smaller than 
the radiation wavelength in the medium. 
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^porous 

Fig. 5. Size dependence of the extinction efficiency fac- 
tors for pseudospheres with inclusions of the same single 
size after averaging of three different targets. Each par- 
ticle contains an equal volume fraction (33.33%) of ACl, 
astrosil and vacuum. The effect of variations of the size of 
inclusions is illustrated. 



The difference between compact and porous particles is 
clearly seen, but the results for porous particles with in- 
clusions and layers look rather similar (excluding, perhaps, 
the height of the first maximum). This behaviour is con- 
firmed by Fig. IHl where other efficiencies, albedo and the 
asymmetry parameter are presented. The largest devia- 
tions occur for the scattering efficiency Qsca (in the range 
a;porous ~ 3 — 10) and the asymmetry parameter g for the 
size parameter a;porous ~ 10. The latter does not seem to 
be an artifact related to small number of angles used in 
our calculations during the averaging over scattering. 

Thus, we can conclude that if materials are "well 
mixed" inside a particle of intermediate porosity, its inter- 
nal structure in form of layers, Rayleigh or non-Rayleigh 
inclusions hardly can be inferred from the transmitted ra- 
diation. In contrast, there is a clear difference between the 
optical properties of compact and porous grains. 

3.2. Very porous particles 

It is commonly accepted that the fraction of vacuum 
in interstellar dust grains can be large. For example, 
very porous particles are often used to model cometary 
grains and dust in protoplanetary discs. Greenberg & 
Hage H1991|l claim that the porosity of dust aggregates 
in comets can be in the range 0.93 < V < 0.98. Their 
conclusion is based on the model of porous grains devel- 
oped by Hage & Greenberg H1990|l who used a volume 
integral equation method similar to the DDA for light 
scattering calculations. A verification of this method had 
been made only for small compact spheres, but the method 
was applied to large and very porous particles. A quali- 
tative agreement between the results obtained with this 
method and from Mie-Garnett calculations was found. 
In both cases the absorption cross section Cabs increased 
and albedo A decreased when the porosity grew. Although 
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Fig. 4. Size dependence of the scattering and absorption efficiency factors, albedo and asymmetry parameter for 
multi- layered spheres. The parameters of particles are the same as in Fig. |3| 
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Fig. 6. Size dependence of the extinction efficiency factors 
for pseudospheres with a size distribution of inclusions. 
Each particle contains an equal volume fraction (33.33%) 
of ACl, astrosil and vacuum. The volume fractions of in- 
clusions of different sizes are approximately the same. 



the validity of these conclusions for particles beyond the 
Rayleigh domain remains unclear, the results of Hage 
& Greenberg (T!J9(T) are frequently used for estimates of 
grain properties in comets (see, e.g.. Mason et al. I2()()l|l . 



Here we analyze in detail highly porous particles in the 
case oi V — 0.9. The results are presented in Figs. 1(^111 
in a manner similar to that used in Figs.|H-|Z| We should 
note the problem connected with the construction of tar- 
gets when the size of inclusions was large. In this case, 
the total number of dipoles from ACl and astrosil was 
~ 36000 — 39000 and they were located in ^ 50 inclu- 
sions which might not always touch others. Nevertheless, 
the changes of the general behaviour of extinction with an 
increase of the number of dipoles in inclusions are clearly 
seen in Fig. |51 The extinction efficiency factors Qoxt be- 
come larger for small values of x and smaller for large ones. 
This essentially deviates from what was observed for par- 
ticles of intermediate porosity (cf. Fig. 01. The growth of 
porosity leads to the disappearance of the first maximum. 
But the curves for particles with inclusions of large sizes 
do not approach the limiting value Qext = 2 defined by 
the "extinction paradox" . This fact should be related to 
special topology of very porous particles with large inclu- 
sions (values of ridip)- Possibly, some inclusions intercept 
a part of light scattered from other inclusions and scatter 
it in the forward direction. This decreases the extinction. 

However, if particles are composed of subparticles of 
different size, their extinction is similar to the usual one 
and the factors Qcxt ^ 2 if a; ^ oo (see Fig.EU- The same 
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Fig. 8. Size dependence of the scattering and absorption efficiency factors, albedo and the asymmetry parameter for 
pseudospheres with inclusions. The parameters of particles are the same as in Fig. 



conclusion is valid for extinction produced by an ensem- 
ble of particles with inclusions of single sizes presented in 
Fig. ^2 This Figure together with Fig. ^] demonstrates 
a very important result: the optical properties of very 
porous layered particles and particles with inclusions are 
similar. Note that in both cases the models of particles 
were constructed in such a manner that the materials in- 
side them were "well mixed" , i.e. the location of inclusions 
in the form of layers or islands is not distinguished. 

This leads to the interesting conclusion that a very sim- 
ple computational model of multi-layered particles seems 
to be of possible use in treating the optics of composite 
grains. 

3.3. Particles of different porosity 

The role of porosity in dust optics can be properly ana- 
lyzed using the normalized cross sections 

^(n) ^ C(porous grain) ^ 

C (compact grain of same mass) 
j>y2/3 g (porous grain) 

Q (compact grain of same mass) 

The quantity C'"' shows how porosity increases or de- 
creases the cross section. Such an investigation has been 



performed by Kriigel & Siebenmorgen H1994|l for absorp- 
tion cross sections of small particles with x < 1. They cal- 
culated the effective optical constants of porous particles 
using, in particular, the Bruggeman mixing rule, and ap- 
plied the Mie theory to get Qabs- Kriigel & Siebenmorgen 
find that the cross sections C^JJ^ increase with V until 
V ^ 0.6 and then decrease. 

Figure El shows the normalized extinction cross sec- 
tions computed for spheres of different porosity. The re- 
sults are plotted in the scale related to the size parameter 
a^compact calculated according to Eq. |j3Jl. Since the extinc- 
tion factors generally decrease when V increases (cf. the 
lower panel of Fig. |31 and Fig. , the values of C^^t are 
greater than unity if the size parameter is smaller than 
~ 1 or larger than ^ 3. Thus, the porosity increases the 
extinction of small and large^ particles. An opposite case 
is observed only in a restricted range of the size param- 
eters Xcompact ~ 1 — 3 whcrc the extinction by compact 
spheres has a maximum (see, e.g.. Fig. Cj). 

As follows from Fig. 1141 such a behaviour of the nor- 
malized extinction cross sections is accompanied by sim- 

^ For very large particles, the normalized cross sections ap- 
proach to asymptotic values C'"' ^ (1 - P)"^/^ (see Eq. @) 
which are equal to 4.64 and 1.59 if V = 0.9 and 0.5, respec- 
tively. 
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Fig. 12. Size dependence of the scattering and absorption efficiency factors, albedo and the asymmetry parameter for 
pseudospheres with inclusions. The parameters of particles are the same as in Fig. 1111 



ilar changes of the scattering and absorption cross sec- 
tions. At the same time, the scattering and absorption ef- 
ficiencies sharply and slightly grow with x for very porous 
grains of large sizes. Note also that both A (beginning 
with a^compact ^2 — 3) and g (for particles of all sizes) in- 
crease with porosity. The behaviour of and A found 
by us is more complicated than predicted by Hage & 
Greenberg (1990). Namely, the growth of porosity leads 
to an increase of C'"' and a decrease of A for very small 
size parameters, and to an increase of both quantities for 
large values of x. There exists also a small interval of inter- 
mediate size parameters where both C^"^ and A decrease. 

Therefore, we can expect larger extinction, scatter- 
ing and absorption of radiation by porous particles with 
radius Ts^compact at wavelengths A ^ 2/37rrs,compact and 
A ^ 27rrs, compact in comparison with compact particles 
of the same mass. At the intermediate wavelengths, the 
compact particles absorb and scatter more radiation. For 
example, the "importance" of compact grains in the pro- 
duction of extinction is larger at the near-UV/visual range 
of wavelengths (0.21/im ^ A ^ 0.63/im) and at the near-IR 
wavelengths (2.1//m ^ A ^ 6.3/^m) if compact = 0.1 /im 
and Ts compact = 1/im, respectively. Note that the latter 
estimates are rather approximate because of the wave- 



length dependence of the refractive index of materials (see 
Voshchinnikov et al. 120041 for more details). 

4. Comparison with Effective Medium Theory 

The EMT is an approach to treat inhomogeneous scat- 
terers as homogeneous particles having an average (effec- 
tive) refractive index. The attempts to find the "best" 
EMT mixing rule continue up to now (e.g., Maron & 
Maron I^UUHl . The EMT is weU described in the recent 
reviews of Sihvola fl9991, Chylek et al. f2000') and papers 
of Spanier & Herman ( 2000) and Kolokolova & Gustaf- 
son H2()()l|l . There are many EMT rules, but besides a few 
ones they are rather similar in principle. Here we give for- 
mulas of the most often used EMT rules for n-component 
mixtures: the Garnett 1)1904(1 and Bruggeman H1935|l rules. 
In the first case, the mixing rule averages the dielectric 
permittivities of inclusion materials et and a "matrix" 
(host) material Sm*-^'' 



/ 



^ m 

e,. + 2e„ 



(5) 



2e„ 



The dielectric permittivity is related to the refractive index 

, ^2 
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Fig. 14. Size dependence of the scattering and absorption normalized cross sections for multi-layered porous spheres 
and pseudospheres with a size distribution of inclusions. The parameters of particles are the same as in Fig. El 




Fig. 15. Size dependence of the efficiency factors (left panel) and their relative errors (right panel) calculated with the 
exact theory for multi-layered spheres and with the Mie theory using four different EMT rules. Multi-layered particles 
contain an equal volume fraction of amorphous carbon (ACl), astrosil and vacuum. The cyclic order of the 18 layers 
is indicated. 



where fi = V^/Vtotai is the volume fraction of the ith ma- 
terial and ScB is the effective permittivity. The expression 
for the Bruggeman H1935|l rule is 



Ecff 



2ec 



= 0. 



(6) 



As an example of a more sophisticated rule, we use the 
"layered-sphere EMT" introduced by Voshchinnikov & 
Mathis (1999). In this case, the effective optical constants 
Eoff are defined as follows (see also Farafonov l2UUUfl : 



Eoff 



l + 2a/V 
1-a/V 



A2 



(7) 



where a is the complex electric polarizability and the co- 
efficients Ai and A2 are obtained as a result of multipli- 
cation of matrices depending on the optical constants and 
volume fractions of layers 



Ai 
A2 



1 



1/3 

-2/(3e„) 




-2/(9/. 




1/3(2^ 



(8) 



The absolute bounds to £eff were given by Wiener 



£off, 



and 



^eff, min 



(9) 



(10) 



These expressions were exactly derived for non-absorbing 
materials but it seems they can be applied to slightly absorbing 
materials too. 
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Fig. 7. The average extinction efficiencies for particles 
with single size inclusions (from Fig. |SJ) and particles with 
size distribution of inclusions (from Fig. |SJ). Each parti- 
cle contains an equal volume fraction (33.33%) of ACl, 
astrosil and vacuum. The thick solid line corresponds to 
compact spheres consisting of ACl and astrosil. For a 
given value of the size parameter, the compact and porous 
particles have the same mass. The thick dashed line shows 
the extinction for layered spheres after averaging over four 
samples presented at the bottom panel of Fig. 13 
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Fig. 9. Size dependence of the extinction efficiency factors 
for very porous pseudospheres with inclusions of different 
size. Each particle contains volume fractions of ACl and 
astrosil equal to about 5%, the porosity V — 0.9. The 
particles are similar to those presented in Fig.|31but have 
larger porosity. 



For any dielectric particle composition and structure, ScS 
cannot lie beyond these limits as long as the microstruc- 
ture dimensions remain small compared with the radia- 
tion wavelength. Note that in the composite grain model 
of Mathis & Whiffcn H1989II the effective refractive index 
was calculated from Eq. @, i.e. a maximum of possible 
refractive indices was taken. 

The general condition of EMT applicability is that the 
size of "inclusions" (in the EMT the particle inhomogene- 
ity is considered in the form of uniformly distributed small 
inclusions) is small in comparison to the wavelength of in- 



size distribution of inclusions 
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Fig. 10. Size dependence of the extinction efficiency fac- 
tors for very porous pseudospheres with a size distribution 
of inclusions. Each particle contains volume fractions of 
ACl and astrosil equal to about 5%, the porosity V = 0.9. 
The particles are similar to those presented in Fig. but 
have larger porosity. 
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Fig. 11. Averaged extinction efficiencies for particles with 
single size inclusions (from Fig. and particles with a 
size distribution of inclusions (from Fig. llO|l . Each particle 
contains volume fractions of ACl and astrosil equal to 
about 5%, the porosity V = 0.9. The particles are similar 
to those presented in Fig.|Z|but have larger porosity. The 
thick dashed line shows the extinction for layered spheres. 



cident radiation (Chylek et al. 12000(1 . The real range of 
applicability of different rules was shown to be nearly the 
same (see, e.g.. Table 4 in Voshchinnikov '2004|) . 

Now let us discuss how different EMT rules can repro- 
duce the optical properties of layered spheres. Figure [T51 
(left panel) shows the extinction efficiency factors com- 
puted with the exact theory for layered spheres and with 
the Mie theory using Garnett^, Bruggeman, and layered- 
sphere mixing rules of the EMT (i.e. with the effective 
refractive indices are equal to rUcs — 1.496 -I- 0.060i, 
TOcff 1-541 + O.OSli and meg = 1.529 + G.080«, respec- 
tively). Wiener's maximum bound is — 1. 604-1-0. 105i. 
Figure El (right panel) demonstrates the relative errors 



Vacuum was adopted as the matrix material. 
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Fig. 13. The normaUzed extinction cross sections for 
multi-layered spheres of different porosity. Open circles 
show the normalized extinction for pseudospheres with a 
size distribution of inclusions (see Fig.lTUIl. 
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Fig. 16. Size dependence of the extinction efficiency fac- 
tors calculated for multi-layered spheres, pseudospheres 
with inclusions and with the Mie theory using three dif- 
ferent EMT rules. The porosity of particles is "P = 0.9. 



of these EMTs. It can be seen that the errors of the 
Bruggeman and layered-sphere rules are of several per- 
cent or better in the considered range of particle sizes. 
The same is generally true for other efficiency factors 
and albedo. An exception is the region after the first 
maximum of the scattering efficiency factor and albedo 



(soporous ~ 6 — 8) where the relative errors may reach up 
to 20%. The largest errors occur for the asymmetry pa- 
rameter, especially for small size parameters. The high ac- 
curacy of the layered-sphere rule in the case of very small 
particles sizes is explained by the fact that it is based on 
the Rayleigh approximation. 

Other rules of the EMT behave like the Garnett and 
Bruggeman rules. The general condition of the EMT appli- 
cability is not fulfilled for layered particles as "inclusions" 
(layers) are not small in comparison with the wavelength 
of incident radiation. However, most rules of the EMT can 
reproduce the optical properties of layered spheres of any 
size, if the number of layers is larger than 15 - 20 (and 



the particles with well mixed inclusions as it was shown 
in Sect. 13.1(1 . This conclusion, however, can be affected by 
the porosity of particles. 

FigureElillustrates the applicability of different EMT 
rules to particles of very high porosity. The cases of other 
efficiency factors, albedo and asymmetry parameter are 
similar. The Figure demonstrates that the layered-sphere 
rule rather well reproduces the optical properties of lay- 
ered spheres as well as the particles with inclusions of 
different sizes (the errors are smaller than 10 - 20% if 
a;porous ^ 15 and V — 0.9). Note that for such particles 
other rules provide acceptable approximations for inter- 
mediate porosity {V ^ 0.5). The Garnett and Bruggeman 
rules together with the Mie theory rather well approxi- 
mate the light scattering by particles with inclusions of 
small sizes. Therefore, all previous models based on the 
standard EMT-Mie calculations are related to particles 
composed of subparticles of very small sizes. If the size 
of subparticles is not small, only the layered-sphere rule 
can be used for the description of the optical properties of 
very porous scatterers. 

5. Concluding remarks 

We consider different (including new) computational ap- 
proaches to calculating the optical properties of composite 
and porous grains that can be used for the interpretation 
of observations of interstellar, circumstellar and cometary 
dust. In our models the particles are represented by multi- 
layered spheres or pseudospheres with inclusions of one or 
different sizes. If the number of layers is small, our model 
of layered spheres coincides with older models of the grains 
having several coatings. For a large (^ 15 — 20) number 
of layers, the model of layered spheres can approximate 
heterogeneous particles consisting of inclusions of differ- 
ent sizes. This gives us a handy way to treat compos- 
ite grains employing a very simple computational model 
of multi-layered particles instead of time-consuming DDA 
calculations. 

We make a careful examination of the optical proper- 
ties of particles of various porosity. Previously, this task 
was solved using the Mie theory for homogeneous spheres 
and effective refractive indices derived from different mix- 
ing rules of the Effective Medium Theory. It is demon- 
strated that this approach gives relatively accurate re- 
sults only if particles have small (Rayleigh) inclusions. 
Otherwise, the approach becomes unacceptable when the 
porosity exceeds ~0.5. An exception is provided by a so- 
phisticated layered-sphere mixing rule, recently suggested 
by Voshchinnikov & Mathis H1999|l . that gives results of 
acceptable accuracy for particles with non- Rayleigh inclu- 
sions as well. Note, however, that our consideration was 
restricted by spheres, non very absorbing materials and 
the integral scattering characteristics but not the differ- 
ential cross sections or elements of the scattering matrix. 

Some astrophysical applications of the model of layered 
grains (in particular, the possibility to reduce the model 
dust-phase abundances) will be presented in a subsequent 
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paper (Voshchinnikov et al. l2(JU4fl . Further development of 
the model of multi-layered particles will involve considera- 
tion of non-spherical inhomogeneous grains (see Farafonov 
et al. I2UU3I for a review). 
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